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THE EVALUATION OF AUTONOMIC RESPONSES: 
TOWARD A GENERAL SOLUTION* 


By 
John I. Lacey 


The Fels Research Institute, 
Antioch College, Yellow Springs, Ohio 


Responses of the central and peripheral structures of the autonomic 
nervous system are ubiquitous. One cannot stimulate the organism, how- 
ever innocuously, without producing some evidence of disturbance of 
autonomic equilibrium. This fact is still not thoroughly integrated with 
our thinking about psychological matters, but the realization is growing 
that restricting consideration of the autonomic nervous system to a 
chapter on ‘‘emotions’’ in our texts and reviews is an underemphasis, or 
at least a misemphasis. 

A brief enumeration of some of the relationships found or suggested 
may serve to illustrate the scope and importance of the measurement of 
autonomic nervous system functions. The activity of the autonomic nerv- 
ous system is related to the adequacy of sensory, sensorimotor, and 
perceptual behavior (E. Calloway and S. V. Thompson, 1953; N. Kleitman, 
S.Titelbaum, and P. Feiveson, 1938; A. Trehub, 1954; J.J. Van Biervliet, 
1894) and to the adequacy of the performance of skilled behavior (T. A. 
Hussman, 1955). Measurements of autonomic responses to the innocuous 
stress of mild exercise (the Schneider test) are related to psychiatric 
status (R. A. McFarland and J. H. Huddleson, 1936), to ‘‘frustration 
tolerance’? (G. B. Strother and D. M. Cook, 1953) and to individual dif- 
ferences in the relative frequency of recall of completed and incompleted 
tasks (T. M. Abel, 1938). Autonomic responses to injections of epineph- 
rine and of mecholyl (the Funkenstein test) are related to recovery rate 
under both psychotherapy and somatotherapies (E. Gellhorn, 1953; L. 
Alexander, 1955). Resting autonomic measurements are related to combat 
performance (M. A. Wenger, 1948), Autonomic responses are being used to 
evaluate the effects of psychotherapy (W. N. Thetford, 1952; O. H. Mowrer, 
1953) and to study the course of the therapeutic interview (A. Di Mascio, 
R. Boyd, M. Greenblatt, and H. C. Solomon, 1955; R. Boyd and A. Di 
Mascio, 1954). These responses are being used in an attempt to under- 
stand the effects of the new tranquilizers such as chlorpromazine (R. A. 
Schneider, 1955). As ‘‘the voice of the unconscious”’ and the main medi- 
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ating structure in the production of psychosomatic dysfunctions, the | 


° . eos | 2 
autonomic nervous system plays a central role in psychosomatic medicine. — 


Suggestively, the autonomic nervous system appears prominently in con-§ 
siderations of ‘‘energy mobilization’’ (E. Duffy, 1951) and the “energetics — 
of behavior’? (G. L. Freeman, 1948); and its activity has been viewed as — 
an important part of both positive and negative feed-back circuits to the 4 
cortical cells, thus controlling, modulating, and terminating cortical | 
activity (C. W. Darrow, 1947). In these last relations, the autonomic ~ 
nervous system loses its traditional role of serving as an “‘objective 
indicator’’ of ‘‘subjective affect,’? and assumes a far more important func- 
tion as a self-regulatory mechanism that is part of the organic government . 
of behavior. : 

This is a broad scope indeed. In an area of such importance, careful © 
and even anxious consideration of the problem of quantifying autonomic . 
responses is desirable. Exploration of this problem immediately raises a 
basic question: In what sense are we justified in saying that individual 
A, or group A, is more or less reactive than individual B, or group B? Or, 
that individual A is more reactive to stimulus a than to stimulus 6? The © 
answer to this question involves at least 3 important issues: the ‘‘law of 
initial values’? and both intraindividual and interstressor differences in 
the organization of autonomic responses. 


eo 


—_" 


The ‘‘Law of Initial Values’’ 


The fundamental problem is to decide on the mathematics to be used. 
That there is a real problem here is evidenced by a large and discomfort- 
ing literature on the measurement of changes in skin resistance. There 
are unsuspected pitfalls awaiting the investigator who simply and ex- 
pediently quantifies skin-resistance responses as an algebraic or per- 
centage change. One of the main sources of difficulty, of course, is that 
the magnitude of the response is related to the prestimulation level. In 
general, it has been found that high autonomic excitation preceding 
stimulation is correlated with low autonomic reactivity upon stimulation. 
In comparing individuals or groups, then, one does not know whether to 
attribute differences in obtained reactivity to pre-existing differences in 
the background level of autonomic excitation, or to reactivity per se. 
Much statistical ingenuity and labor has ‘een expended in demonstrating 
this phenomenon and in suggesting and testing mathematical transforma- 
tions of the raw data of skin-resistance response with the hope of arriv- 
ing at a meaningful and base-free measure of skin resistance change that 
would be generally applicable (J. P. Seward and G. H. Seward, 1935; C. 
W. Darrow, 1937; O. L. Lacey and P. S, Siegel, 1949; E. A. Haggard, 
1945, 1949; A. S. Paintal, 1951; D. N. Elliott and E. G. Singer, 1953; M. 
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A. Wenger and O. C. Irwin, 1936). A bewildering array of transformations 
has resulted: the logarithm of the change in conductance; the change in 
the logarithm of conductance; logarithm of the change in resistance; the 
ratio of the logarithm of the change plus a constant divided by the rest- 
ing level; the percentage change in resistance; the square root of con- 
ductance; the percentage any given ohmic decrease is of the maximum 
ohmic decrease obtained, and others. 

It is significant to note that the transformations have not been 
rooted in physiological fact or theory but have emerged entirely as a 
tesult of statistical manipulations. Moreover, they have turned out to be 
Specific to the original data on which they were developed, and investi- 
gators who actually go to the trouble of testing them with newly obtained 
data seem likely to find that they do not work (A. S. Paintal, 1951; P. H. 
Venables, 1955; R. A. Champion, 1950; R. A. Terry, 1953). 

Strangely enough, intensive consideration has been given to the 
problem only with regard to skin resistance. The same problem, however, 
appears in the measurement of heart rate, blood pressure, blood flow, 
muscle potentials, blood-sugar levels, and so on. A new and a more 
general solution is needed. 

I begin my approach to a more general solution with the bald con- 
tention that a basic physiological error is involved in the common prac- 
tice of quantifying the reactions of autonomically innervated structures 
by computing algebraic or percentage changes, and that additional error 
is committed by seeking arbitrary mathematical transformations that ir- 
rationally (by which I mean empirically, and without regard to physiologi- 
cal theory) remove the dependence of such computed changes upon pre- 
stimulus levels of physiological functioning in a specific experimental 
Situation. I propose that such procedures in effect ignore the basic fact 
that one of the main functions of the autonomic nervous system and as- 
sociated endocrine glands is to maintain a homeostatic norm. I further 
propose that the properties of an adequate statistical model for the eval- 
uation of autonomic responses must correspond to the real properties of 
the universe of events we wish to study. In what follows, I shall first 
attempt to delineate these ‘‘real properties,’’ which affect our measuring 
operations, and I shall then attempt to develop an adequate statistical 
model. 

Homeostatic restraint of response. A huge number of experimental 
observations suggest that any induced excitation or inhibition of an 
autonomically innervated structure instantly initiates a series of changes 
that serve to nullify the disturbance. There is an important corollary 
which somehow has escaped being set forth in italics in discussions of 
»sychophysiological response mechanisms. The corollary is this: the 
recorded autonomic response 1s a ee both of the induced magnitude 
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of autonomic activation (as it would be seen in the absence of contrary 
changes) and of the promptness and vigor of secondarily induced autonom=— 
ic changes that serve to restrain and limit the effects of the initial “/ 
turbance. Let us be bold about it. This corollary should be made part of 
the implicit background of every psychophysiological west 
conclusions should be drawn in such investigations without explicitl i] 
taking the corollary into account. . 
Darrow (C. W. Darrow, 1943), in a brilliant review that should be a 
classic reference for psychophysiologists, attacks the problem from a_ 
different but revealing viewpoint. ‘‘The problem,’’ he states, ‘‘is liter-_ 
ally to determine the weight on either side of a ‘balance’ when that of 
neither side is known.’? Darrow considers ways of circumventing the 
“mutually antagonistic action of the two branches of the autonomic 
nervous system’’ by surgical, pharmacological, and experimental-sta- 
tistical techniques. He desires to find means of making valid inferences 
concerning the individual roles of the sympathetic and parasympathetic 
systems in the production of a peripheral response of a given magnitude. 
We desire to find means of evaluating the final resultant reaction, in the 
belief that it is the final resultant that is of importance to a variety of 
problems central in psychophysiology and psychosomatic medicine. | 


The 2 problems, however, are intimately related. Evidence is avail-— 
able that suggests that it is precisely because of certain characteristics 
of the ‘‘mutually antagonistic action’’ of the sympathetic and para- 
sympathetic that we should expect a relationship between prestimulus 


physiological level and stimulus-produced physiological changes. 

Homeostatic restraint of cardiovascular response. We can marshal 
this evidence best by an elementary review of cardiovascular dynamics, | 
emphasizing those aspects that clarify the problem of quantifying a pe- 
tipheral response. While the relevant material is available in many books 
and articles, a most recent authoritative presentation has been made by 
Heymans (C. Heymans, 1950). 

As is well known, peripheral resistance and cardiac output are con- 
tinuously and with great sensitivity modified to maintain arterial blood 
pressure within relatively narrow limits. An important part of the mecha- 
nism underlying this control lies in the aortic arch and in the carotid 
Sinus, innervated by the so-called buffer nerves (carotid sinus nerve and 
cardio-aortic depressor nerve). As the arterial pressure rises, the pressure- 
Sensitive receptors in these strategically located structures are stimu- 
lated and, through their connection with the vagus, may initiate wide- 
Spread changes: a decrease in the rate of impulse formation in the sino- 
auricular node; a diminution of the force of contraction of the auricles 
and ventricles; a reduction of the duration of systole; a decrease in the 
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conductivity of the impulses through the different parts of the heart; and 
peripheral vasodilation. These changes obviously reduce cardiac output 
and peripheral resistance, and thus lower the blood pressure. With this 
compensatory lowering of blood pressure, the inverse changes occur, 
with the result that cardiac output and peripheral resistance increase, 
with a consequent rise in arterial pressure. This continuous feed-back 
‘mechanism maintains the homeostatic norm of blood pressure and heart 
tate, 

The precision and sensitivity of the mechanism are strikingly ex- 
emplified by 3 observations: (1) in the carotid sinus of the dog, a change 
of pressure of as little as 1 mm. Hg is sufficient to initiate compensatory 
changes in blood pressure; (2) at each systole there occurs an outburst 
of action potentials over the carotid sinus nerves and, with each diastole, 
a cessation of activity, at low to moderate arterial pressures; while at 
high pressures the impulses form a continuous train (C. Heymans, 1950; 
D. W. Bronk, and G., Stella, 1932; D. W. Bronk, 1935); (3) not only is the 
absolute level of intrasinus pressure an adequate stimulus for the pres- 
sure receptors, but so is the rate of change of that pressure (D. W. Bronk 
and G, Stella, 1932). 

Of utmost importance to the argument of this paper are the following 
observations by Samaan (A. Samaan, 1935). 


(1) While the reflex regulation of cardiac output is maintained both 
by the vagus and the cardiac sympathetic accelerator nerves, changes in 
vagal tonus are far more important and significant. Thus, cardiac slowing 
(dogs and cats) produced by moderate vagal stimulation masked the effect 
of simultaneous strong cardioaccelerating stimulation. For example, to 
cite a typical result, if there were simultaneously applied to the pe- 
tipheral cut ends of the appropriate nerves a sympathetic stimulus capable 
of increasing heart rate by 41 beats/minute and a vagal stimulus capable 
of decreasing the heart rate by 34 beats/minute, the net result is a de- 
crease of 30 beats/minute. 

(2) Sympathetic excitation potentiates vagal inhibition. Thus, simul- 
taneous stimulation of the cardioaccelerator nerves or small doses of 
epinephrine make the heart even more sensitive to vagal stimulation. That 
epinephrine sensitizes carotid sinus mechanisms has been abundantly 
verified (C. W. Darrow, 1943). 

(3) The vagal effect is far more prompt than the acceleratory effect. 
Cardiac slowing upon vagal stimulation appears at the first beat that 
follows stimulation. Cardiac acceleration, however, does not make its 
appearance until a minimum of 2% seconds after the onset of stimulation, 
and the maximal acceleratory effect does not appear until a minimum of 
7 seconds has elapsed. —s 
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7 
Putting all these effects together-the continuous and sensitiveld 
reflex regulation; the greater importance of parasympathetic than of 
sympathetic control; and, most important, both the potentiation of vagal 
influence by simultaneous sympathetic excitation and the greater rapidity 
of vagal effect-we are clearly led to expect that the magnitude of the — 
absolute or relative change that can occur in an intact organism will be 
negatively related to the physiological level at the onset of cs 
Put in other terms, the more powerful influence of parasympathetic in- 
hibitor mechanisms, and their potentiation by sympathetic excitation, — 
suggest the following interpretation of the dynamics of a cai ae 
response. As the prestimulus level of functioning increases, there is a 
disproportionately greater homeostatic restraint, both in increased magni-_ 
tude and decreased latency and, as the magnitude of induced activation — 
increases, there is a disproportionately greater increment of counter-_ 
reaction. : 
Homeostatic restraint in other organ systems. Such effects probably 
are not limited to the cardiovascular system. The secretion of epinephrine 
is self-limiting; the presence of increasing amounts of circulating epi- 
nephrine in the blood increasingly inhibits the secretion of further epi- — 
nephrine and inhibits transmission across ganglionic synapses (A. Mar- ; 
razzi, 1939). In a similar manner, circulating epinephrine limits the re-— 
activity of other sympathetic functions (C. W. Darrow, 1943), including . 
sweat-gland reactivity. The carotid sinus mechanism itself extends con- — 
trol over a wide variety of somatic functions. Literature reviewed by | 
Darrow (C. W. Darrow, 1943) shows that increased intracarotid pressure © 
blocks respiration, halts shivering, reduces or eliminates skeletal re- 
flexes, and decreases muscle action potentials. | 
Many details need to be worked out, indeed, but enough is at hand to — 
suggest that the negative relationship of stimulus-induced physiological — 
change to prestimulus physiological level is not simply an arithmetical 
artifact of a restricted range of response left available to the organism 
whose prestimulus functioning is high. It 1s a phenomenon of physio- — 
logical homeostasis and must be dealt with as such, If statistical trans- 
formations are sought in order to derive base-free measures of response, 


they must not be arbitrary ones. They must be related rationally to the 
physiological phenomenon. 


Experimental observations in intact organisms. The theory outlined 
above is not completely unchallenged by experimental data. That homeo- 
Static restraint of response may fail to be effective on occasion is only 
to be expected. Very recently, Brown and his collaborators have reported 
experiments on anesthetized dogs that show little or no effect of carotid 
Sinus reflexes in moderating systolic blood-pressure response to injected 
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pinephrine, norepinephrine, and histamine (G.C. Boxill and R.V. Brown, 
953; R. V. Brown and J. G, Hilton, 1954; J. G. Hilton and R. V. Brown, 
954, 1955). It becomes important, then, for the psychophysiologist to 
est the generality of occurrence of the phenomena of homeostatic re- 
traint in intact human organisms. 

The first approach to a general recognition of the relation of auto- 
omic changes to prestimulation level was made by Wilder (J. Wilder, 
950) in a paper that has appeared only, I believe, as an abstract. Wilder 
ave the phenomenon a useful name: the ‘‘law of initial values.’’ Work- 
ing, apparently, with a variety of pharmacologic agents in serial experi- 
ments with man, Wilder found in ‘‘75 to 85 per cent of the experiments’’ 
that the magnitude of increase of physiological changes decreases as the 
initial level increases. Indeed, at high initial levels he often found no 
increase and ‘‘quite often’’ a ‘‘paradoxical’’ decrease. These are sig- 
nificant observations, and we have verified them repeatedly in our own 
laboratory. A technique of quantifying autonomic responses certainly 
should not obscure these phenomena. 

Wilder also noted another significant fact that, with inhibitory agents, 
the opposite rule held. The higher the initial level of autonomic excita- 
tion, the greater the decrease in the physiological function under ex- 
amination. 

How well does this ‘‘law’’ of initial values hold up? In TABLE 1 are 
presented the results of several of our studies involving a variety of 
populations, stresses, and physiological variables. The details of the 
experiments are published elsewhere (J. I. Lacey and R. Van Lehn, 1952; 
a I. Lacey, D.E.Bateman, and R. Van Lehn, 1953). In TABLE 1 the 
children, subjects in the Fels longitudinal study of growth and develop- 
ment, range in age from 6 to 18, and are of both sexes. The mothers, also 
subjects in the longitudinal study, range in age from 23 to 57. The col- 
lege males are a much more homogeneous group, being freshmen at 
Antioch College, Yellow Springs, Ohio. Our technique involves first a 
15-minute period of relaxation to stabilize physiological levels. The sub- 
ject is then warned that in 1 minute he will have to put his foot into a 
cold-water bath, or will have to do mental arithmetic, or hyperventilation, 
ot letter association. This results in an increment of autonomic activity 
which we call the alerting reaction. After 1 minute of such anticipation, 
the subject is then required to undergo the indicated stress, Another 
increment of activity results that we call the stress reaction. The sub- 
ject is then allowed to recover to his previous resting level before the 
next alerting period begins, and so on seriatim. For college males, there 
was no alerting period. 

One characteristic of our treatment of data should be noted. Each 
physiological variable differs, of course, in its mode of response. These 
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variables have different latencies. One may appear as a smooth mono- 
phasic response, another as a polyphasic response, etc. To put all 
variables on an equal footing, and to get clearly interpretable data, we 
measure the maximum sympatheticlike activity. Thus, for heart-rate alert- 
ing reaction, we measure and average the 6 fastest heart cycles in the 
Jast minute of rest, and then the 6 fastest heart cycles in the minute of 
anticipation of stress. For heart-rate stress reaction, the latter figure is 
Compared with the 6 fastest heart rates during stress. 

I shall shortly attack the use of correlation coefficients in this 
Situation, and shall show that they give an inadequate, and occasionally 
an erroneous, picture of the true state of affairs. With this proviso, one 
Can expediently get some impression of the generality of the law of 
initial values. Considerable sampling variation is evident in TABLE l, 
especially and surprisingly for systolic blood pressure (but compare the 
findings of Brown et al., cited above). There is no room for doubt, how- 
ever, that initial level is negatively related to percentage change over 
a wide variety of stresses and populations. The same relationship holds 
for algebraic change, with one startling exception. For palmar conduc- 
tance (measured in log micro-ohms), we secure significant positive cor- 
relations. One possible reason for this will become apparent in the dis- 
cussion below concerning the complexities of meaning of correlation co- 
efficients between initial level and a measure of change. Before leaving 
this table, it is well to note that these are linear correlations, and that 
bivariate normality is not assured. In many instances (see below), curvi- 
linear correlations would reveal relations of increased magnitude. 

In FicuREs 1 and 2 are shown the results of serial experiments with 
college males. These are subjects participating in an hour-long session 
devoted to the study of semantic conditioning and generalization of un- 
conscious anxiety (J. I. Lacey and R. L. Smith, 1954; J. I. Lacey, ewes 
Smith, and A. Green, 1955). Heart-rate responses repeatedly are elicited 
at 45 second intervals by combined motor activity (tapping a telegraph 
key) and verbal activity (chain association to stimulus words), For 28 of 
such episodes of reaction spread through the hour, the coefficient of 
correlation between initial value (6 fastest cardiac cycles during the last 
15 seconds of a 30-second intertrial rest period) and change to stress 
level (6 fastest cardiac cycles during 15 seconds of verbal-motor activity) 
was computed for each of 21 subjects. The results again strongly support 


the law of initial values. 


The Use of Regression Models 


What can we do about this dependence of change upon initial level? 
Desirable and rational experimental arrangements are obvious. In 
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serial experiments with man, impose each succeeding stressful crisoay 
upon the same initial level. This is possible, but many times ex rena 
difficult to attain, In group experiments, allow each ag ear to relax ort 
recover to a predetermined level, constant for all subjects, before | 
ministering the stress. This is manifestly impossible and, in itself, it is | 
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CORRELATION BETWEEN PRESTIMULUS VALUE 
AND ABSOLUTE CHANGE IN HEART RATE 


FIGURE 1. Frequency distribution for 21 cases of the individual product- 
moment coefficients of correlation between prestimulus heart rates and the 
algebraic increments of heart rate due to combined verbal-motor activity. The 
cardiac activity during each episode of stimulation is measured by averaging 
the 6 fastest cardiac rates in 15 seconds of combined activity. From this value 
is subtracted the average of the 6 fastest cardiac rates in the last 15 seconds 
of the preceding 30-second intertrial rest period, to yield the algebraic increment 
in ‘‘maximum sympatheticlike activity.’? Each correlation is based on 28 serially 


elicited reactions. Note that all correlations are negative, and that the median 
correlation is —0,73. : 
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CORRELATION BETWEEN PRESTIMULUS VALUE 
AND PERCENTAGE CHANGE IN HEART RATE 


FIGURE 2. Frequency distribution for 21 cases of the individual product- 
oment coefficients of correlation between prestimulus heart rates and the per- 
tage increments of heart rate due to combined verbal-motor activity. Pre- 
timulus and stimulated levels of heart rate are determined as in FIGURE 1. 
ach correlation is again based on 28 serially elicited reactions. Note that all 
strelations are negative, and that the median correlation is —0.75. 


“technique that poses interpretative difficulties. The ideal solution 
ould be to have extensive statistical norms for varying and well-defined 
spulations and stresses. These norms would present separate frequency 
istributions of changes for varying initial levels. A given individual’s 
Sponse, then, could be referred to the frequency distribution of re- 
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sponses obtained at an initial value identical with that shown by the 
individual being studied, and his response could be expressed as the 4 
number of standard deviations away from the mean response. This ideal, 1 
of course, will never be attained with any broad generality. Given a well- 
defined experimental situation, however, and a not too prohibitively large 
sample of a well-defined population, we can attain a satisfactory and 
useful approximation to this ideal solution by utilizing statistical tec 
niques of so-called regression analysis. A best-fitting curve, in the 
least-squares sense, relating the x-arrays (initial levels) to the means o } 
the y-arrays (responses), becomes the locus of estimated average per- 

centage or average algebraic changes at varying initial levels, and the 
standard deviation about this ‘‘line of regression’’ is an estimate of th ‘ 
standard deviation for any distribution of responses at any initial level. 


If the data behave in certain lawful and simple ways, and if certain 
restrictive assumptions are met, sensitive and fairly precise inferences” 
can be drawn, using only a knowledge of the means and standard 7 
tions of the initial levels and the responses of our sample, and the cor- 
relation between initial level and response. What are the requirements 
(beyond adequate sampling of a defined population)? 

First, the best-fitting curve should be linear. If it is not, and- 
curvilinear relationship fits the data significantly better than a linear 
one, we should be faced with the arduous and difficult task of decidin 
exactly what the form of the equation is and the equally arduous task of 
deciding on the parameters and constants of the equation. If we attempt 
to avoid this by first grouping our subjects into arrays of varying initial 
levels and then computing means and standard deviations for each such 
subgroup of individuals, we shall be losing considerable sensitivity and 
ease of application for the following reasons. 

For the linear case, we have the simple and well-known result that 
Yz =TxyXz, where y, is the y-value on the regression line that corre- 
sponds to a given x, value (as the z subscript implies, x and y values 
are measured in units of sigma deviations for the entire sample), and 
r,, is the product-moment correlation between initial level and response. 
That is, knowing an individual’s initial level, we can solve the equation 
to estimate the theoretical average response at that initial level. This’ 
regression equation can be derived from ungrouped data, and we can makel 
our inferences sensitive within one-integer class-intervals of initial 
level. That is, we can estimate mean response for individuals with initial 
levels of heart rate of 80 beats/minute and a different mean response for 
individuals with initial levels of 81 beats/minute, and so on. The only 
limitation lies in the sensitivity and precision of the original data. In 
the curvilinear case, however, no such Simple result is likely to occur. 
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e shall not, in general, have a simple or rationally acceptable regression 
quation to solve; because of this we shall have to resort to coarse 
touping of initial levels, and we shall have to compute means and 
tandard deviations of the responses for each gross class-interval of 
nitial level, The lack of precision and reliability of such a procedure is 
bvious. 

Second, the curve should be homoscedastic. That is, the standard 
eviation of responses should be the same from x-array to x-array, vary- 
ing only at random and within the limits of sampling error. If the curve is 
oth linear and homoscedastic, we again have a simple and well-known 
a the standard deviation of responses, at any initial level, is given 

oy /1-r2, , where %y is the standard deviation of the total distribu- 
a of y wae refinement of securing a completely unbiased estimate by 
multiplying °y ./1-r2 by \ (W-1)/(W-2) is unnecessary for reasonably 
large N’s). 

Third, the y variable should be normally distributed over all x-arrays 
and within each x-array. 

These conditions of linearity of regression, homoscedasticity, and 
normality, then, define the situation where the product-moment coefficient 
of correlation is maximally useful. 

In the course of these investigations I tested the assumptions of 
linearity, homoscedasticity, and normality in 306 separate plots involving 
all the data in rasBLe 1 and additional data. I consider the results of 
only historical interest because, as I shall show shortly, regression 
models involving percentage or absolute change are inappropriate, and 
they contain hidden traps. Nevertheless, they are worth summarizing. 


| (1) The normality assumption is frequently violated. Not only are 
distributions of absolute change or relative change not normal, but fre- 
quently no empirical transformation can be found that will normalize the 
distributions. Sixty per cent of the distributions of percentage change 
were nonnormally distributed. Approximately 40 per cent of the distribu- 
‘ions of algebraic change were nonnormal. 

(2) Moreover, the regression of response upon the initial level is 
often decisively nonlinear. Approximately 16 per cent of the plots show 
significant curvilinearity. 

(3) The assumption of equal variance in the x-arrays was acceptable 
in every plot where normality and linearity were also satisfied. 

On these grounds alone we should reject the applicability of regres- 
sion models involving measurements of change, but there are even more 
convincing reasons for abandoning the approach. Interestingly enough, a 
sonsideration of these reasons directly suggests the appropriate and 
yorkable approach. 

ae 
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Regression techniques involving percentage or algebraic change in. j 
volve (1) redundancy, and (2) ‘‘spurious’’ correlations. . 


Redundancy. The element of redundancy is obvious. Consider a sub > 
group of individuals at a given initial level, for example, at 110 mm. Hg } 
blood pressure. Subject them to a stressor agent and record the maximum 
blood pressure reading. Since all individuals have started with the same 
initial level, we get no more information by computing percentage or 
algebraic increment from this initial level as a base than we get by — 


gests that the most direct and meaningful technique is to remove the re- © 
gression of stress level upon base level. No intermediate calculation of 
change is required. | 

‘‘Spuriousness.’’ Let x be the initial level and y be the stress level. | 
Then A = (y-x) is the algebraic increment of autonomic activity, and the 
ratio @ = (y-x)/x is the proportional increment. In correlating x with (y-x)) _ 
or x with (y-x)/x, the element x is common to both variables being cor- _ 
related. Some intriguing and dangerous complications follow this introduc- } 
tion of ‘‘spurious’? communality. As shown in the appendix to this paper, _ 
written by Edmund Churchill, the correlations between initial level and” 
either percentage or algebraic change are directly a function of the cor- 7. 
relation between initial level and stress level, and of the variances of 
those levels. . 


For algebraic change, an exact identity is given by 
, %y? (1 ~r2,) | 
rmrzl- Eman ac, meas where the symbols have their usual 
y -2r,y, yx+x q 

meaning, 


Note some consequences, If %? = %?, and the basic correlation 
between x and y is 0, then r,) will be equal to —,707! In other words, 
even if there is no correlation of initial level with stress level, we shall 
get some oorrelation of initial level with algebraic change, depending 
upon the relationship of the standard deviations of the 2 levels. As 
the ratio of the component standard deviations varies, so will the cor- 
relations. If the ratio °x/ y is greater than r,y, the correlation of x with 
(y-x) will be negative; given the same value for txy, the correlation of x . 
with (y-x) will be positive if %&/%y is less than tx» The resultant base- © 
change correlations are then extremely difficult to interpret, being a. 
result simultaneously of the correlation between the component measures © 
of level, and of the relationship between their standard deviations. The, 
Situation is far too complex to allow a simple and meaningful physiologi- 
cal interpretation, me 
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_ For percentage change, an approximate identity is given by 


+1 71+3V? 
‘= —————_—_————__ ,,_ where x, y, and @ are defined as 


1+8V7 1-2 


_ HO 
1+3v’? [v 2 
ws” 


bove, V is the coefficient of variation of x, and W is the coefficient of 
atiation of y. Here the complexities are of the same kind, except that a 
ominant role is played by the coefficient of variation, rather than the 
standard deviation, r,g taking on various values as a function of the re- 
lationship between the ratio of the 2 coefficients of variation, and the 
basic correlation r,,, that is, as a function of [(V/W) —r]. Both the sign 
and the value of the correlation between x and 0 depend on the relative 
value of r,, and the ratio of the 2 coefficients of variation. Some of the 
misleading consequences of this situation are detailed in the appendix. 

Let us summarize. The basic data are the initial level and the stress 
level. Computing percentage or absolute change is then a redundant step. 
Somputing the correlation between such changes and the initial level 
ntroduces a maze of undesirable complexities. 

We are then left with the simple recourse of dealing directly with 
he invariably high positive relationship (found for the group, not, as we 
shall see below, for the individual) between initial level and stress level. 
*or the same data where so much difficulty was found in establishing 
satisfactory conditions for the use of the regression model using per- 
centage or absolute change, I found no obstacle in this simpler and 
nore direct application. All distributions of both initial and stress level 
yere normal, by statistical test, either in their raw form, or by using the 
yiologically familiar and defensible logarithmic transformation (J. H. 
zaddum, 1945). All bivariate distributions yielded linear regressions and 
lomescedasticity by appropriate statistical tests. 

Our proposed general measure of autonomic reaction, then, is simply 
ittained, by repetitive solution, ‘for each individual, of the equation: 


Vz — Xz buy 


Tene)” 


\utonomic Lability Score =50 + 10 , where x, and y, are 


he individual’s initial and stress levels, respectively, expressed in 
inits of sigma deviation in the total sample, t,, is the correlation for 
he sample between initial level and stress level, and the constants 10 
ind 50 translate the resulting scores to a distribution with a mean of 50 
nd a standard deviation of 10. This last step is simply an arithmetic 
‘convenience in that it dispenses with negative numbers and provides a 


; ne 
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system of units that makes it easy to compare an individual’s reactio 
in any number of physiological variables simultaneously. 

The interpretation of the final result is clear. Consider an individual] 
whose autonomic lability score for blood-pressure response to the col ie 
pressor test is 50. The meaning is that the level which blood pressute | 
reached during stress was exactly that level predicted from his initia 
level. If his response were 60, the meaning would be that the level his: 
blood pressure reached during stress was 1 standard deviation (theo- 
retically estimated) above the mean level (theoretically estimated) at-- 
tained by those of his peers who had the same initial level. 


Autonomic Lability Scores: Some Results and Problems 


Certain important characteristics of the scores should be summ 
ized: 

First, they are, of course, free of dependence on base level, and are 
normally distributed if the assumptions underlying regression analysi 


are fulfilled (as they have been to date in our studies). . 
Second, in accordance with physiological theory, the regression line 


autonomic activation is induced upon a prestimulus background of high 
excitation. 

Third, equal absolute or percentage changes (if we compute them), 
Starting from different initial levels, do not receive equal autonomi¢ 


a given change gets larger and larger autonomic lability scores as the 
initial level increases, the amount by which the autonomic lability 
Scores increase being determined by the correlation that reveals the 
degree of homeostatic restraint for the given variable, population, anc 
Stress, These points are illustrated in TABLES 2 and 3, using data fo 
heart-rate responses obtained serially in man. 

Fourth, at zero correlation between initial level and stress level 
(a result we have sometimes obtained in serial experiments), the auto- 
nomic lability score is identical, except for the translation to a dis 
tribution with a mean of 50 and a sigma of 10, with the sigma score for 
Stress level. This, we feel, is as it should be. If there is no homeostatic 
relation to take into account, then the level of activation under stress i 
the only important variable. As the correlation between initial level and 
stress level increases, the autonomic lability score is increasingly 
divergent from the corresponding sigma score for stress level. 

In general, autonomic lability scores correlate well with absolute or 
percentage change per se. In one study (J. I. Lacey, D. E. Bateman, and 
R. Van Lehn, 1953), we obtained correlations ranging from 0.74 to 0.93 
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ith a median value in the high 80’s. These are high correlations, and 
one might wonder if the game were worth the candle. One example will 
Serve to show how misleading absolute or percentage changes can be, 
and how illuminating autonomic lability scores can be. 

In one unpublished study, we investigated maturational changes in 
patterns of autonomic activation by the cold-pressor test, using both al- 
gebraic and percentage change, for a sample of 64 children ranging in 
age from 6 to 18 years. We established some tidy results. For example, 
younger children tended to overreact in diastolic blood pressure and to 
underreact in heart rate, whereas older children showed the reverse 
pattern. 

The correlation with age of a pattern-index score (defined in the 
teference given immediately above), indicating the degree to which 
diastolic blood-pressure reactivity exceeded heart-rate reactivity, was 
-0.60, using a percentage measure of change and —0.43, using an al- 
gebraic measure of change. Confidence levels were obviously very satis- 
factory. Our confidence, however, was shaken when we found in our own 
data the marked age changes in resting level that other investigators 
have found. As children grow older, their diastolic blood pressure in- 
creases and their heart rate decreases. As children grow older, then, we 
might expect (on the intriguing hypothesis that longitudinal experiments 
over a period of years will reveal, in the individual, the same result seen 
in short-time serial experiments) that it would be harder to evoke a 
marked increase in blood pressure, but easier to evoke a marked increase 
in heart rate. When we used autonomic lability scores, the alleged age 
change in the pattern of autonomic reactivity disappeared completely, 
and the correlation with age dropped to 0.07. The age change in pre- 
stimulus level, then, completely accounted for the apparent change in 
reactivity. 

It is apparent, therefore, that although the correlation may be high 
between absolute or percentage change on the one hand and autonomic 
lability scores on the other, enough variance in the data remains to be 
attributed to base level to mislead us in interpretation. Our results to 
date indicate that whenever autonomic lability scores show a significant 
variation with our independent variable (age, sex, intensity of shock, 
reinforcement schedule, Rorschach responses, and the like), absolute and 
percentage changes do also, The autonomic lability scores, however, 
have always yielded clearer trends and greater interpretive clarity. On 
the other hand, absolute or percentage changes may yield results that 
do not hold up when the effect of base level is removed, as in the ex- 
ee just given. ee 
‘ Intraindividual autonomic lability scores. We attach special interest 


ee 
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and importance to the technique of regression analysis applied to serial 
experiments with the single individual, and we have successfully applie 
intraindividual autonomic lability scores in studies of semantic con= 
- ditioning (J. I. Lacey and R. L. Smith, 1954; J. I. Lacey, R. L. snl 
and A. Green, 1955). 

Certain special problems, however, arise in this application. Our 
experiments, in the citations already given, involved only 28 serially 
elicited responses. This yields an inadequate number of degrees of free- 
dom to enable proper evaluation of normality, linearity, and homoscedas- | 
ticity. Inspection of each case, however, seemed to reveal acceptable 
properties, and we simply went ahead with the analysis. In current ex-_ 
periments on the same problem we have 52 serially elicited responses, 
and we shall soon be in a position to make appropriate statistical tests. : 
In the meantime other problems have arisen, and one intriguing and 
significant by-product has appeared. 

The problem is: frequency distributions of initial levels and of 
stress levels for individuals sometimes are markedly nonnormal. For 
these individuals, the usually trustworthy logarithmic transformation — 
failed, but a square-root transformation worked. For others, no reason-_ 
able or simple transformation would work. Such individuals remain a 
puzzle to us. To bypass the difficulty, we have resorted to ‘‘artificial”? — 
normalization by means of McCall’s T-scores (P. O. Johnson, 1949), 
While not overly pleased with this somewhat irrational maneuver, the 
results were clear and reasonable in that the final experimental results 
for individuals so treated do not differ from those secured from individu- 
als whose distributions are transformable to the normal form. + 

This procedure, however, recommended itself as an expedient one to | 
use for all individuals, since it would remove the necessity of laborious- 
ly testing for normality of distribution in each individual case. No in- 
formation seemed to exist, however, on the accuracy of such ‘artificial’? | 
normalization, especially as a function of the size of the sample. In a 
previous study (J. I. Lacey, D. E. Bateman and R. Van Lehn, 1953), we 
showed that the autonomic lability scores using the McCall transforma- 
tion correlated from 0.91 to 0.99 with the scores resulting from the use 
of logarithmically normalized data, for an N of 85. Since then we have 
undertaken an empirical study, using hypothetical and actual data for 
heart rate, blood flow, and palmar conductance variables for which it was 
known that logarithmic or square-root transformations were effective. The 
nature of the distributions may be seen in TaBLeEs 4 and5. The degrees | 
of skewness vary over a wide range, and of kurtosis over a narrow range, 
characteristic of the type of data with which we deal. As can be seen 
from the tables, the raw distributions were such that the ratio &,/SEg, 
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TABLE 4 


DEGREES OF SKEWNESS (G,) AND KURTOSIS (G,) FOR HYPOTHETICAL AND 
ACTUAL DISTRIBUTIONS USED TO COMPARE MCCALL’S T-SCORING 
METHOD WITH. A NORMALIZING SQUARE-ROOT TRANSFORMATION* 


Untransformed distribution Transformed distribution 


*For convenience of presentation, the figures are rounded to 1 decimal place. 


enables rejection of the hypotheses that the distributions are normal, 
_ whereas the transformed distributions are acceptably normal (P. O. John- 
‘son, 1949). From these distributions 2 sets of scores were constructed. 
One set of scores, in T-score form (with a mean of 50 and a standard 
deviation of 10), was constructed using the mathematical transformation 
that was known to normalize the raw data. The second set of scores was 
emade up ‘‘artificially,’’ using the McCall technique. The results of com- 
_ paring these 2 sets of scores are shown in FicurREs 3, 4, 5, and 6. It 
_ may be seen that ‘‘artificially’’ normalized scores converge rapidly upon 


; TABLE 5 

-. 

"DEGREES OF SKEWNESS (G,) AND KURTOSIS (G,;) FOR HYPOTHETICAL AND 
ACTUAL DISTRIBUTIONS USED TO COMPARE MCCALL’S T-SCORING 


METHOD WITH A NORMALIZING LOGARITHMIC TRANSFORMATION* 


Untransformed distribution’ Transformed distribution 


Pays? 


*For convenience of presentation, the sieures are rounded to 1 decimal place where 
possible, . 
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SQUARE-ROOT TRANSFORMATION 


CORRELATION 


0 10 20 30 40 50 60 70 80 90 100 110 IZ09 


NUMBER OF CASES 


™“ 


FIGURE 3. Coefficient of correlation, as a function of sample size, between 
physiological levels expressed in O units based on a square-root transformation 
that has been demonstrated (see TABLE 4) to normalize the distribution, and 
levels expressed in 0 units based on McCall’s transformation. A curve has been 
drawn by visual inspection to suggest the form of the relationship. Note that the 
results of the 2 transformations rapidly converge on each other. 


scores normalized by a mathematical transformation. N’s of 40 give al- 
most as good results as N’s of 120. 
These figures, however, do not make one point clear. The greatest 
differences between the 2 sets of scores invariably occur at the ex- 
tremes of the distributions. The reason is that McCall’s technique 
truncates the distributions at both ends. The results as a whole suggest | 
that the use of McCall’s T-scores can be safely recommended for data 
such as we obtain, even with relatively small N’s, although one has to 
be careful in certain applications in evaluating extremely high OF @x- 
tremely low reactions. z 
It should be noted, before leaving this topic, that iitraindividuall 

it 
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regression analysis of data obtained in serial experiments requires that 
several responses be obtained at each of many prestimulus values and 
that account be taken of the correlation of the components of response 
with time. In many lengthy serial experiments, it will be found that a 
process of adaptation sets in that is revealed partially by a monotonically 
(or almost monotonically) increasing or decreasing prestimulus level. 
Thus, afterwhile, each succeeding response in the series is elicited from 
Successively higher or lower prestimulus values. This phenomenon of 
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FIGURE 4. Average difference, as a function of sample size, between 
ysiological levels expressed in 0 units based on a square-root transformation 
at has been demonstrated (TABLE 4) to normalize the distribution, and levels 
xpressed in 0 units based on McCall’s transformation. A curve has been drawn 
visual inspection to suggest the form of the relationship. 
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FIGURE 5S. Coefficient of correlation, as a function of sample size, be- 
tween physiological levels expressed in 0 units based on a logarithmic transform 
mation that has been demonstrated (TABLE 5) to normalize the distribution, and 
levels expressed in 0 units based on McCall’s transformation. A curve has been 
drawn by visual inspection to suggest the form of the relationship. Note that the 
results of the 2 transformations rapidly converge on each other. 


almost perfect correlation between prestimulus value and time is seen 
frequently in conditioning experiments when skin-resistance responses 
are measured. The technique of intraindividual regression analysis is 
inapplicable in such circumstances, of course. Some other technique of 
analysis will have to be devised to suit each experiment. | 
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Individual differences in homeostatic efficiency. It may have oc- 
‘curred to the reader that a physiological interpretation of the ‘‘law’’ of 
initial values, such as we have proposed here, generates a new question 
‘and a new vantage point from which individual differences in autonomic 
“teactivity may be surveyed. The question was foreshadowed in FIGURES 
“1 and 2, where tremendous variation among individuals can be seen in 
the magnitude of the negative correlation between prestimulus values of 
heart rate and stimulus-produced changes. Does this variation reflect in- 
dividual differences in homeostatic efficiency, or is it mere sampling 
variation? 
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FIGURE 6. Average difference, as a function of sample size, between 
physiological levels expressed in 0 units based on a logarithmic transformation 
nat has been demonstrated (TABLE 5) to normalize the distribution, and levels 
expressed in 0 units based on McCall’s transformation. A curve has been drawn 
by visual inspection to suggest the form of the relationship. 
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A fresh grasp on this question is provided by. gos tesults obteiia 
in as-yet-unpublished studies on semantic conditioning, which are cm 
tinuations of the studies already cited. In these new experiments well 
elicit serially 52 heart-rate responses by combined mental-motor activity 
in an anxiety-conditioning experiment. FicurReE 7 presents a freque . 
distribution, for 56 subjects, of individually computed correlations 


OF CASES 
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CORRELATION BETWEEN PRESTIMULUS & 
STIMULUS LEVELS OF HEART RATE 


| 
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FIGURE 7. Frequency distribution for 56 cases of the individual product- 
moment coefficients of correlation between prestimulus values of heart rate and 
stimulated values of heart rate due to combined verbal-motor activity. Distribu- 


tions of heart rates for each individual have been normalized by McCall’s trans- 
formation. Each correlation is based on 52 serially elicited reactions, The 
cardiac activity during each episode of stimulation is measured by averaging the 
6 fastest cardiac rates in 15 seconds of combined activity, and is correlated 
with the average of the 6 fastest beats in the last 10 seconds of the preceding 
30-second intertrial rest period. Note the large individual differences in home- 


ostatic efficiency. The rank-order split-half reliability of the correlations is 
+0.60, uncorrected (see text). 
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etween initial level and stress level using McCall’s normalizing tech- 
ique. Again, tremendous variation is shown, with the correlations falling 
nto class intervals whose midpoints vary from —0.05 to +0.85. Some in- 
ividuals seem like automatons. The level their heart rate reaches during 
xperimental stimulation is closely related to the prestimulus heart rate. 
her individuals seem much more variable; their cardiac rates during 
timulation are unrelated to the background level of autonomic excita- 
ion, each succeeding episode of verbal-motor activity is responded to 
almost uniquely. 

To study the matter further we selected a sample of 20 individuals, 
all of*whom were in the same experimental group (there are 4 such 
soups, with systematic variations in experimental treatment among them). 
[he individually computed correlations of these individuals ranged from 
+0.19 to +0.79. The pool of 52 serial reactions for each individual was 
split into 2 halves, which were equated for temporal position and seman- 
ic connotation. Two correlations were then computed for each subject, 
| for each equivalent subgroup of 26 responses. These pairs of cor- 
elations were then themselves correlated by the rank-difference tech- 
lique. The resulting correlation was +0.60, significant at beyond the 
).01 level. Within this single experimental situation, then, the variation 
f homeostatic correlations (as we may call them) is not entirely random. 
Phere are reliable and systematic individual differences. 

- In these data, there is a possibility of some circularity of reasoning. 
Phe individual correlations were obtained in a study of the growth and 
lecay of conditioned autonomic responses. Perhaps it is precisely those 
ndividuals who exhibited the greatest conditioning and generalization 
vho show low correlations. This would follow if the increasing autonomic 
tivation seen as the conditioning proceeded was very marked in absolute 
erms. If, for example, before conditioning, heart-rate changes were of 
he order of magnitude of 5 beats/minute, but during and after condition- 
ng were of the order of 30 beats/minute, a low correlation between pre- 
timulus and stimulus values might be produced as a consequence of the 
ery psychological process we are studying. That this potential circular- 
ty of reasoning does not exist is suggested by the fact that the quanti- 
ative variation of response is not nearly so marked as in the hypotheti- 
al example cited above, and by the details of the arithmetical pro- 
edure by which we secure indices of conditioning and generalization 
see references given). We attacked the problem directly, however, by 
ialyzing curves of conditioning and generalization as a function of the 
omeostatic correlation.’’ No statistically significant differences were 
‘ained. Whatever trend there was, moreover, was in a direction that 
ssured us there is no circularity in our data. It was the individuals 


a 
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with higher correlations who seemed to condition and generalize note 
promptly and extensively. 4 
Further analysis showed that these differences in homeostatic f. 
ficiency were not related to (1) the variation among individuals in “7 
all level of heart rates, or (2) the variation among individuals in th 
variability of prestimulus levels of heart rate. Thus, the magnitude ° 
the homeostatic correlation is not to be attributed to variations in the 
background level of autonomic excitation (such that high correlatio 
would be found in individuals whose heart rates were high, thus mor 
urgently invoking homeostatic restraint). ; 
We tentatively conclude that the apparent and large individual dif- 
ferences in homeostatic efficiency are not attributable to chance varia 
tion. They seem real, and it seems likely that the set of individual dif 
ferences so revealed will turn out to be of some importance. The source 
of the variation and the physiological and behavioral consequences ne 
to be ascertained. 


All this, taken as a whole, gives a first answer to our origina 
question. We are justified in comparing the reactivity of groups, or 
individuals, or of a given individual to varying stimuli, only to the de 
gree to which we realize and account for the fact that the autonomic re 
action we record is a complex manifestation of opposing influences. 
tational solution to one important problem posed by the existence 
such homeostatic mechanisms — the law of initial values — may be a 
proached by using the regression of stress level of functioning upon p 
stimulus level of functioning. 

This is not the whole story, however. A second and most provoca 
tive complication is produced by the fact, only recently established ex 
perimentally, that autonomic responses are systematically organized int 
patterns of activation. Limitations of space forbid any but the briefes 
review of the facts in this area of investigation. 


The Patterning of Autonomic Responses 

Interindividual response specificity. The correlations among aaa 
nomic responses to stress are usually rather low, indicating either a lack 
of communality or of reliability. In 5 separate experiments, only 3 of 
which have been published as yet (J. I. Lacey, 1950; J. I. Lacey and 
R. Van Lehn, 1952; J. I. Lacey, D.E. Bateman, and R. Van Lehn, 1953), 
we have found strong evidence that the lack of correlation is due be 
systematic and idiosyncratic patterns of response. The main facts are 
well illustrated in ricure 8, which shows the responses of a 45-year- 
old woman to 3 different stresses. These are autonomic lability scores 


” 
= 4 
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franged in T-score form with a mean of 50 and a standard deviation of 
0. The graph illustrates 2 phenomena: 

- (1) A marked intraindividual differentiation of autonomic response. 
hether we call this person markedly overreactive, markedly under- 
active, or of average reactivity, depends on the function being measured. 

(2) This individual shows what we have called relative response 
pecificity. Whatever the stress, she tends to respond with the same 
ierarchy of activation. 

Statistical and experimental analyses of various sorts show in 
arious ways that this is a reliable phenomenon. Individuals tend to re- 
spond to different stresses, which impose varying physiological and 
»sychological demands on the organism, with reproducible patterns of 
autonomic activation. Over a period of 4 years, moreover, the patterns 
9f response to the cold-pressor test are highly reproducible. What is not 
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65 
tw 
ow 60 
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< Legend: 
>> Cold Pressor 
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30 —-—- Letter Ass'n. 
HEART- HEART SYSTOLIC PULSE PALMAR DIASTOLIC DIGITAL 
RATE RATE BLOOD PRESSURE CONDUCTANCE BLOOD BLOOD 
VARIABILITY PRESSURE PRESSURE FLOW 


PHYSIOLOGICAL VARIABLE 


_ FIGURE 8. Profile of the autonomic activation for a 45-year-old woman 
hho was subjected serially to 3 stresses: cold-pressor test, mental arithmetic, 
d letter association. Autonomic reactions are expressed by autonomic lability 
cores (see text). Note the marked intraindividual differentiation of autonomic 
sponse. The woman is markedly overreactive in heart rate and heart-rate 
ariability, for example, but markedly underreactive in digital blood flow. Note 
Iso that she exhibits ‘‘relativeresponse specificity.’? The general pattern of 
ctivation is the same from stress to stress. 
Be 
7 
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yet clear is whether these results are a statistical expression of a | 


a 
‘ 


tendency common to all people to yield idiosyncratic patterns of re 
sponse, or whether there are quantitative variations in this tenden v3 
towards differentiation and stereotypy of patterns of activation. Pt 
general biological grounds, and on the basis of our data to date, we fa I 
the latter hypothesis, and we are undertaking experiments to test if: 
However this latter hypothesis tums out, it is clear that intraindividual 
differentiation of autonomic response is a factor to be considered. Both} 
the direction and degree of differentiation among individuals depends or 
the physiological function measured. Using one physiological measut 4 
we should claim that A is less reactive than B; using another, we should: 
claim the reverse. Corroboration of this principle is provided by Te 7 
(R. A. Terry, 1953) and, for the case of individuals with frank psycho-4 
somatic disorders, by Malmo and his collaborators (R. B. Malmo, C.. 
Shagass, and F, H. Davis, 1950), who independently formulated a closely ' 
related principle of ‘‘symptom specificity.’”’ | 
One possible way to rescue a concept of over-all autonomic reactive > 
ity is indicated in a previous publication (J. I. Lacey, D. E. Bateman, | 
and R. Van Lehn, 1952). We were there able to show a significant rela 
tionship between the Rorschach form-color index of ‘‘emotionality’’ an 
autonomic responses to stress, only if the individual’s reactivity was; 
assayed in terms of his own pattern of response. Two individuals, in: 
other words, who exhibited a maximum autonomic lability score of 62) 
let us say, were called equally reactive, although 1 individual showed 
his maximal score of 62 in heart rate, and the other in palmar conduct | 
ance. This specific relationship, however, was not reproducible in an 
unpublished study by Vaughn Crandall and myself on a different popula: | 
tion of individuals. Obviously, much study and thought will have to be | 
given to the problem of evaluating such patterns. 
Interstressor response specificity. The pattern of autonomic activa j 
tion is different in different affective states. Thus, Ax (A. Ax, 1953) 
has been able to show physiological differentiation between responses : 
he labeled ‘‘anger’’ and ‘‘fear.’? Funkenstein and his collaborators have 
aoe asieencee in the patterns of arousal in responses labeled 
anger-out,’’ ‘‘anger-in,’? and ‘‘anxiety’’ (D. H. Funkenstein, S. H. King 
and M. Drolette, 1954), The relationship between these findings and the 
principles of ‘‘symptom specificity’? and ‘‘relative response specificity’ | 
is still unclear. A direct and quantitative attack on the problem of what 
proportion of the variation in response is to be attributed to individue 
differences in the organization of the autonomic nervous system and what 
Ptoportion is to be attributed to the nature of the individual’s psycho 
Sain mechanisms of coping with specific stressor episodes has yet to 
e made, s 


a 
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Summary and Conclusions 


Now let us retrace, in a summary, the long and rather involved path 
own which we have come. We have shown, for a variety of populations 
nd stresses and for cardiovascular and sudomotor variables, that the 
agnitude of stimulus-produced activation of the autonomic nervous 
ystem, in general, is related negatively to the prestimulus level of physio- 
ical function. We have argued that this is not simply an arithmetical 
ttifact of there being but a limited range of response left available to an 
rganism already functioning at high levels of autonomic excitation, but 
hat it is a physiological consequence of, and deducible from, the facts 
f homeostatic restraint of response. We believe that the same mecha- 
sms operate for variables other than those of the cardiovascular and 
judomotor, and we have briefly touched upon the relevant literature. 

We have proposed that arbitrary and empirical statistical transforma- 
ions causing this phenomenon to disappear within a given body of ex- 
erimental data are irrational procedures, and that the rational procedure 
s to establish statistical norms, that is, frequency distributions of 
esponses at varying prestimulus levels of functioning for varying popu- 
ations and stresses. Since this is manifestly impossible, we propose 
hat the statistical best approximation to this ideal solution is attained 
y use of the technique of regression analysis. 

We then attempted to demonstrate that removing the regression of 
Igebraic or percentage response on initial level was spurious, redundant, 
nd awkward. We proposed then to analyze only the regression of stress 
svel on initial level. The scores resulting from such a regression analy- 
is we called autonomic lability scores. 

This technique proved adequate in all respects, and it yielded clear- 
‘*para- 


y interpretable results. The technique, moreover, freely allows 
oxical’’ phenomena to be exhibited, that is, an apparent lack of auto- 
omic increment or the appearance of autonomic decrement when experi- 
‘ental stimulation is imposed upon high background autonomic excita- 
ion. Such ‘‘paradoxical’’ events may receive high autonomic lability 
cores in the proposed regression model. 

- We then illustrated the utility of autonomic lability scores in adding 
larity of interpretation to some data on maturational changes in patterns 
f autonomic activation. 

The concept and the measure were then extended to the individual 
ase. One difficulty in this application appeared in that nonnormal dis- 
ibutions (within individuals) were found. The applicability of an “arti- 
cial’? normalizing technique (McCall’s T-scores) was investigated by a 
mited empirical study. The results suggested that such an “‘artificial’’ 
schnique may be used safely with the sort of data we obtain. — 


oe 
“ 
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It was shown that both the viewpoint and the data generated r 
questions and new phenomena: Are there reliable individual differenc es 
in homeostatic efficiency revealed by the coefficient of correlation be 
tween prestimulus and stimulus levels of autonomic excitation? T 
data presented tentatively suggest an affirmative answer. 

Finally, we briefly considered the fact that autonomic activatit 
does not occur as a quantitatively undifferentiated and mass discharg 
Instead, it is patterned. There are both reliable interindividual and re 
liable interstressor differences in the patterns of activation. , 

It is clear, then, that there are 3 aspects to the answer to our origi+ 
nal question. To evaluate autonomic responsivity we must (1) take the 
law of initial values into account, (2) choose autonomic variables th 
are relevant to the individual, and (3) choose variables that are releva 
to the behavioral sequence we are studying. For the first problem, the 
elements of a general concept applicable to the formulation of a general 
solution are at hand. In facing the problems of patterned activation off 
the autonomic nervous system, we are at the very edge of knowledge. A} 
vast area of ignorance faces us, to challenge the efforts of psychologist, 
physiologist, and physician, of clinician and laboratory scientist. 


APPENDIX 
By Edmund Churchill 


Department of Mathematics, Antioch College, Yellow Springs, Ohio 


To find the correlation between the variable x and the variable d 
(y-x), with x and y expressed as raw scores, let x, %, ¥, % and can 
the parameters of the bivariate normal distribution of x and y, and let j 
stand for the ‘‘expected’’ or mean value. Then, from the standard e 
pression for the correlation coefficient: 


Oren E [x(y-x)] — E (x) E(y-x) 


Oy a 
Now E[x(y-x)] = E(xy — x?) = E(xy) — E(x), 
E(x) = x, and 
E(y-x) = E(y) — E(x) = (y-x). 
Then, E[x(y-x)] — E(x) E(y-x) = E(xy) — E(x?) -x y+x- 
= (E(xy) -¥ 7] - (E(x -=7] 
= Le Oy Oy = o2 . 


And, as is well known, 


G, 5 . 
(ex) = NF -2tey Ge oy + OF 
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nce: 
ine Oy Oe = 2. 
tf xy Ux Yy x 
OS a 
@ i 2 
Ce hae ot, (CG. + 0, 


t, letting d = °x/%, and substituting %9 for % in FormULA 2, we get 
d 


i) (3) 2,4 a Fxy Fy — Py 


| F — Ir, 02 cee Pie 


; = poem se Bk." of 
f “Sear ear 
Sa N a 


ince the denominator is “y-x)/7y, and consequently positive, the cor- 
ect value of the square root is the positive one. Hence, if r,,< ¢, that is, 
E xy < %/%, r, will be negative; if 4, >, % will be positive. If 
~ = 9, tA will clearly not be zero, but will equal —¢ / \1 +? or 
@./ \o2 +92, A related formula for the correlation between x and y+ x 
s given by Snedecor (G. Snedecor, 1946, p. 163). Returning now to 
‘ORMULA 2, we have, by a series of simple steps: 


pet eta Parke Rat 
ue 
> 
ll 


on eet ee 


eo ae. 


og — te) 
SEL ya ORR ES RES ee ene which 


2 2 
Oty 036, +0; 


he result given in the preceding article. The sign of the square root 
iven by the sign of the quantity (rx —¢). 


We 
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To discuss the correlation between x and the variable 0 = (y-x)/x3 
we proceed as follows, Let x and y be as above, let V = o0,/x be the cox 
efficient of variation of x, and W = o,/y be the coefficient of variatior 
of y. We assume V and W to be sufficiently small to permit ignoring 
powers higher than the fourth of V, W, and their products, in what fol 
lows. It is probable that with V and W no greater than 25 per cent, the 


resulting error in r,9 will be no more than 2 per cent. Then, as we at 
show, the correlation r, g is given by: , ; 
ae | 
Av 1 + 3V? ; 
(4) 4t.69 = Se 

eres Wes = 

ae pened Set ae ; 

erate ae ~) 
: ; 


where A = -1 if (V —1r,x,W) is positive, +1 if (V — r,,W) is negative. // 
more approximate formula is given by Snedecor (G. Snedecor, 1946, p. 16 4) 
The sign of the middle term of the denominator in Snedecor’s formula 
however, appears to be in error. ] 


Clearly, it is the quantity S — fy} tather than oe itself whiel 


plays the dominant role in determining the value of r,9. This fact can | e 
pointed up by considering a few special cases. | 
For esate, if r,, = 1, r,9 will be pasar gar to i 


most eeNajaly not zero, marie being equal to 


-V1+3V? 
a ~ —________ 
1+8V? | (Ww)? 1+W? 
1+3¥° (Vv)? V2 


When V =W, rx@ becomes 


To simplify the derivation of FORMULA 4, let E represent, as abov 
the expected or  yhcee let x =x (1+ u), y=Vy 7 +v). Then E(u) = 
E(v) = 0, E(u*) =V E(v”) = w?. Further, r- = , because u is 
linear transform of x, baa v of y. 
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_ Since (1 + wu) is also a linear transform of x and since, because: 


0 = (y-x)/x = y/x -1 aA AL ey) -1, 6 and (1 + v)/(1 + u) 
(x) +a) > 


‘transforms of each other, the correlation between x and @ is identical 
th that between (1 + uw) and (1 + v)/(1 +o). 


nce: ; (1 +v) 
(5) z[a+9-G+2 | -E(l+wWE (1+u) 


r.6 = = 
1 +u) (1 + v) 
(1 + u) 


1 but 2 of these terms can be evaluated quickly: 


(1 + v) 


]-#a+9-1 


For the other terms, we make use of the following expansions, ob- 
ined by use of the binominal expansion theorem: 


(9) 1/1+u)=(1+u) =1-u+u -u +u'-... 


10). 1/(1 +0)? =(1+a7 =1--204+3u07 —4u° +5at-... 


_ Thus, using FoRMULA 9 we get: 


an 2023 -2[ ag F)]-e[aroa-aewrweut...) 


=E [d-wru? ui su. .)4+(v-uv+u?v-u3v+utv —-.. | 

= E(1)-E (u) +E (u’) —E(u*) +E (u*) +E (v) -E (uv) +E (uv) - 
E(u5y)+E(u4v)-... 

=1-0+V7-0+3V*+0-s,, VW+0-3r,,V W+0-... 


Siry-43V —r,,.VW—3r,,V W-... 


The evaluation of the individual expected values above (and of those 
ow) follows the formulas given by Kendall (M. G. Kendall, 1943, p,.89, 
tion 3.15; and p. 129, formula Diol). 


= 
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Using rorMULA LO, we get 


(12) £ fae E(l+v) cy =E [savy (1-2u43u?4u45u~.0), 


=E (1-2u4+3u?4u°45u*-.. .) +2(v-2uv+3u?v4u3vi5uty—.. ) 
+ (V -2uv?43u?v?-4u3v745utv7—.. .) 


= E (1) —2E(u) + 3E(u?) — 4E(u*) + 5E(u*) + 2E(v) - 4E(uv) ¥ 
6E(u*¥) - 8E(u?v) + 10E(u*v) + A oe 2E(uv?) + 3E(u?v4 
~4B(a'v") + SEC ee | 
S10 FSV a Ores 3V 405 eae Wt ~ 8 -3r,, VW + 0+: 
W? 0.4 3 +225) VW? ~0 5:3 (1 +4ry) VW? | 


=143V? +154 —4r,, VW — 247, V°W + W? +3(1 +2r%,) VW 
+15 (14+4r2,) VW? + . 


Combining FormuLAs 11 and 12, we find (ignoring all terms fo 
which the combined exponents of V and W exceed 4) that: 


2 
2(l+v)_,f(l+v) (1+ v) 
Te! eS ue cea [ wen 


=143V? +15V4—41, VW —24r,,V°W +W? +3(1 +2r2,) V2W? + ah 
[ 14+2V? +7V* -2r.y VW trey VW? — 8rxy vw] 
=V’ +W? —2r,, VW + 8V* —16r,, V°W+3V7W? +522, V7W? | 
= (1+8V7) (V7—2r,y VW) + (14+3V7)W? + 522, VW? 3 
= (14+8V7) (V?—2r,y VW +12,W?) +(14+3V?) (1—r2, Ww? 
=(1+8V") (V-r,, WY + (1+3V7) (1-12 yw? 


Substituting the values of rormucas 6, 7, 8, 11, and 13 in vom 
LA 5, we get (ignoring higher powers of V and W): 


- 2 te 
nee 1 [i sas3v4 (V?—tyy VW) 


V | 1+8V%) (V=r, WY +(1 +34) (1 ~ 2 yw? 


4 = 


—(1+3V7) (V—r,, W) 


(1+8V7) (V-r,, WY +(1 BV) 0-12, Ww? : 
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iding numerator and denominator by (/1+3V2) (V -r,, W) gives 
(14)’ a = A 1 + 3V2 

2 Rae 
14+8V : oe ee 
1+3V’ (V-2,,W)’ 


mre A=-1 if V >r,W, A=+1 if V < tay W. The function of A is to 
i FORMULA 14’the same sign as FORMULA 14. The positive values 
> always assumed for the square roots. 
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